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Abstract: We consider cosmological parameters estimation in the presence of a non-zero 
isocurvature contribution in the primordial perturbations. A previous analysis showed that 
even a tiny amount of isocurvature perturbation, if not accounted for, could affect standard 
rulers calibration from Cosmic Microwave Background observations such as those provided 
by the Planck mission, affect Baryon Acoustic Oscillations interpretation, and introduce 
biases in the recovered dark energy properties that are larger than forecasted statistical 
errors from future surveys. 

Extending on this work, here we adopt a general fiducial cosmology which includes a 
varying dark energy equation of state parameter and curvature. Beside Baryon Acoustic 
Oscillations measurements, we include the information from the shape of the galaxy power 
spectrum and consider a joint analysis of a Planck-like Cosmic Microwave Background 
probe and a future, space-based, Large Scale Structure probe not too dissimilar from 
recently proposed surveys. 

We find that this allows one to break the degeneracies that affect the Cosmic Microwave 
Background and Baryon Acoustic Oscillations combination. As a result, most of the cos- 
mological parameter systematic biases arising from an incorrect assumption on the isocur- 
vature fraction parameter /j so , become negligible with respect to the statistical errors. We 
find that the Cosmic Microwave Background and Large Scale Structure combination gives 
a statistical error <j{fi so ) ~ 0.008, even when curvature and a varying dark energy equa- 
tion of state are included, which is smaller than the error obtained from Cosmic Microwave 
Background alone when flatness and cosmological constant are assumed. These results con- 
firm the synergy and complementarity between Cosmic Microwave Background and Large 
Scale Structure, and the great potential of future and planned galaxy surveys. 
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1. Introduction 

The standard cosmological model assumes adiabatic initial conditions for the perturbations. 
This simple adiabatic picture is well-motivated as it is predicted by the simplest inflationary 
single- field models [1, 2], and so far it provides an excellent fit to current data (e.g., [3]). 
There is however no a priori reason to discard more general initial conditions involving 
entropy isocurvature perturbations as they can arise from several different mechanisms 
e.g., multi-field inflation, [4, 5, 6, 7, 8, 9, 10], neutrino isocurvature perturbations [11], 
axion dark matter [12, 13, 14, 15, 16, 17, 18, 19, 20] or the curvaton scenario [21, 22, 23]. 
Pure isocurvature models have been observationally excluded [24, 25, 26, 27], but current 
observations still allow for an admixture of adiabatic and isocurvature contributions [3, 28, 
29, 30, 31, 32]. 

To date, Cosmic Microwave Background (CMB) data have been the main source of in- 
formation not just about cosmological parameters but also about the nature of cosmological 
initial conditions. Relaxing the assumption of adiabaticity and allowing an isocurvature 
contribution introduces new degeneracies in the parameters space which weaken consid- 
erably cosmological constraints e.g., [34, 30, 40, 35, 33, 36, 37]. Up to now, the CMB 
temperature power spectrum alone has not been able to break the degeneracy between the 
nature of initial perturbations (i.e. the amount and properties of an isocurvature com- 
ponent) and cosmological parameters; adding external data sets somewhat alleviates this 



issue for some degeneracy directions e.g., [28, 31]. As shown in [38], the precision polariza- 
tion measurements of future and on-going CMB experiments like Planck will be crucial to 
lift such degeneracies. 

Ref. [39] investigated the effect of relaxing the hypothesis of purely adiabatic initial 
conditions and found that even a tiny isocurvature fraction contribution, if not accounted 
for, can lead to an incorrect determination of the cosmological parameters and can affect 
the standard ruler calibration from the CMB. In fact, the presence of an isocurvature 
component changes the shape and the location of the CMB acoustic peaks, mimicking 
the effect of parameters such as Q mo h 2 , Hq and w (see also [30] and references therein). 
Ref. [39] found that this has a crucial effect on "standard rulers" , like the sound horizon at 
radiation drag, inferred from CMB observations. This is of relevance for the next generation 
of galaxy surveys which aims at probing with high accuracy the late time expansion and 
thus the nature of dark energy by means of Baryon Acoustic Oscillations (BAO) at low 
redshift (z < 2). In fact, for Large Scale Structure (LSS), even a tiny isocurvature fraction 
contribution, if not accounted for, can affect the BAO estimation, introducing systematics 
in the BAO observables which can be larger than the statistical errors expected from future 
surveys. 

BAO measurements from galaxy surveys have the potential to be an unprecedented 
powerful probe of the low redshift Universe, and thus of dark energy ([55, 54, 56] and 
references therein). The BAO in the primordial photon-baryon fluid, responsible for the 
characteristic peak structure of the CMB power spectrum, leave, in fact, an imprint in the 
large scale matter distribution and, at each redshift, their physical properties are related 
to the size of the sound horizon r s (zd) at radiation drag (i.e., when baryons were released 
from the photons). Since the CMB in principle can precisely provide the standard ruler, 
i.e. r s (zd), by measuring the BAO location at low redshift with LSS surveys, it is possible 
to probe the expansion history, i.e., the Hubble parameter H(z), and the angular diameter 
distance Da(z) at different redshifts, and thus the dark energy properties. 

It has been found that the effect of a number of theoretical systematics such as non- 
linearities, bias etc., on the determination of the BAO location can be minimized to the 
point that BAO are one of the key observables of the next generation dark energy exper- 
iments. However, it is important to investigate and test the robustness of the method to 
all possible theoretical uncertainties. 

The analysis in Ref. [39] focused on the forecasts from a CMB-only Planck-like exper- 
iment to quantify the impact of the isocurvature contribution on cosmological parameter 
estimation and degeneracies with r s {z ( i) and, therefore, with BAO measurements. Here we 
extend the work to a joint analysis of a CMB Planck-like experiment combined to a future 
space-based LSS survey with characteristics not too dissimilar from the proposed Euclid 
mission 1 . One of the key points is to see whether adding the full cosmological information 
provided by future LSS surveys could eventually break the degeneracies introduced by the 
presence of an isocurvature fraction in the initial conditions. In the present work in fact, 
we exploit, in the LSS case, not only the information enclosed in the BAO positions, but in- 



1 http://sci.esa.int/euclid 



elude the galaxy power spectrum shape adopting the so-called "P(fc)-method marginalised 
over growth information" [41, 65]. Moreover we consider spatial curvature and a more 
general case of dark energy with a varying equation of state parametrized by wo and w a . 
The rest of the paper is organized as follows. In Sec. 2 we review theory and notation 
for isocurvature perturbations. In Sec. 3 we outline our methodology, paying particular 
attention to the method used to compute expected constraints from LSS surveys, including 
not only the BAO location measurements but also the power spectrum shape. We present 
our results in Sec. 4 and the conclusions in Sec. 5. 

2. Isocurvature: theory and notation 

Even if pure isocurvature models have been ruled out [24, 25, 26, 27], current observations 
allow for mixed adiabatic and isocurvature contributions [3, 28, 29, 30, 31, 32]. 

In fact, besides primordial adiabatic perturbations, there can exist the so-called isocur- 
vature or entropy perturbations. These are associated to fluctuations in number density 
between different components of the cosmological plasma in the early Universe, well before 
photon decoupling, and are generated by stress fluctuations through the causal redistribu- 
tion of matter under energy-momentum conservation. 

Density perturbations are then produced by non-adiabatic pressure (entropy) pertur- 
bations (see [42] for a detailed description). 

The entropy perturbation between two particle species (i.e., fluid components) can be 
written in terms of the density contrast 5i = 5pi/pi and the equation of state parameter 
Wi as 2 : 

Sxr = j£- - j£-, (2.1) 

1 + wx 1 + wy 

which quantifies the variation in the particle number densities between two different species, 

and it is equivalent to Sxy = Snx/nx — 5ny jny. 

The most general description of a primordial perturbation accounts for 5 non-decaying 

(regular) modes corresponding to each wavenumber: an adiabatic (AD) growing mode, a 

baryon isocurvature mode (BI), a cold dark matter isocurvature mode (CDI), a neutrino 

density mode (NID) and a neutrino velocity mode (NIV) (see [11] and references therein). 

For the dark matter component 5n c /n c = 5p c /p c , so that the CDM isocurvature mode can 

be written as: 

3 

S C = 5 C - -Sj, (2.2) 

where 5x = Spx/px is the energy density contrast of the X particle species. The baryon 
(b) and the neutrino (v) isocurvature modes take the form: 

S b = S b - -^ (2.3) 

S„ = \d u - ^ 7 . (2.4) 



2 Recall that the continuity equation, which follows from the energy- momentum conservation, takes the 
form p = -3H(P + p) = -3Hp{l + w). 



For the adiabatic mode S c = S u = St, = 0. 

One of the crucial distinctions between the adiabatic and the isocurvature models 
relies on the behavior of the fluctuations at very early time, during horizon crossing at the 
epoch of inflation (or an analogous model for the very early Universe). In the adiabatic 
case, constant density perturbations are present initially and imply a constant curvature 
on super-horizon scales, while in the pure isocurvature case there are not initial density 
fluctuations which, instead, are created from stresses in the radiation-matter component. 

An example of isocurvature model is represented by the curvaton scenario [21, 22, 23], 
which is an alternative to the single field inflationary model. This model relies on the infla- 
ton, the light scalar field that dominates the background density during inflation, and the 
curvaton field a which, decaying after inflation, seeds the observable cosmological pertur- 
bations. The initial conditions then correspond to purely entropy primordial fluctuations, 
namely isocurvature perturbations, because the curvaton field practically does not con- 
tribute to metric perturbations. Besides the fact that current data are compatible with 
this theoretical picture, the curvaton model attracted growing attention because it pre- 
dicts primordial non-gaussianity features of the local type in the spectrum of primordial 
perturbations. 

2.1 Notation 

A common parametrization for the isocurvature perturbations is given by [43]: 

(S 2 \ x l 2 

hS ° ' (Ka d ) 1 ' 2 ' { } 

defined as the ratio between the entropy S and the curvature (adiabatic) 1Z perturba- 
tions evaluated during the radiation epoch and at a pivot scale k$. In our case we set 
0.002 Mpc -1 , as often done in the literature [44, 3]. The correlation coefficient can be then 
defined in terms of an angle, A^ , such that: 

COS A *<> = (V 2 V1/2/C2 U/2 - ( 2 - 6 ) 

V^rad) ' \^rad) ' 

Throughout this paper we will use a sign convention such that fi so < will correspond to 
correlated modes. 

2.2 Isocurvature and CMB 

In general the two-point correlation function or power spectra for the adiabatic mode, the 
isocurvature mode and their cross-correlation can be described by two amplitudes, one 
correlation angle and three independent spectral indices (n a d, n iso , n cor ), so that, in the 
case of the CMB, the respective algular power spectra are: 



^ = /f[<(^(0 



(2.7) 



cr=- l^m so (k)?(^) (2.8) 
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Here @^ d (k) and Ql so (k) are the radiation transfer functions for adiabatic and isocurva- 
ture perturbations that describe how an initial perturbation evolved to a temperature or 
polarization anisotropy multipole I. The total angular power spectrum takes then the form: 

C t = {Tll ad )[Cf + fl Cr + 2f tso cosA k0 Cn- (2.10) 

From Eq. (2.10) it is clear that, for small isocurvature fractions, the main isocurvature 
contribution come from the mixing term coefficient (2 fi so cos A& ). 

2.3 Isocurvature and LSS 

The power spectra of both adiabatic and isocurvature perturbations, as well as their cross- 
correlation, can be parametrized by three power laws with three amplitudes and three 
spectral indices n a d, ^iso and n cor : 



k 3 k 3 / k \ n «d-! 



",.. 



h 3 h 3 / h 

k 3 
A-Ks(k) = 7T~ 2 (7^rad<Srad) 

jl3 / u \ n cor+\{n a d+n iso )-l 

= — %AB cosAz,,, f — ) 

As before, 7^ ra( j stands for the curvature perturbation, and 5 ra( j = (<5 C — 35 7 /4) ra( j for 
the CDI perturbation evaluated during the radiation epoch at a pivot scale k$ such that: 
A = (Kl d )^ and B = (5 r 2 ad ) 1/2 - 

The last of Eqs. (2.11) refers to the extra-correlation that can be generated by the 
partial conversion of isocurvature into adiabatic perturbations after the end of inflation. 

Since we are interested in the effect of a model with a mixed adiabatic+isocurvature 
CDI contribution to the LSS, it is useful to give the explicit formula of the shape of the 
linear matter power spectrum (e.g. [29]): 



P(k) = (A 2 + B 2 )[{\ - a) P ad (k) + a P iso (k) + 2/3^(1 - a) P cor (k)\ . (2.12) 

Here P ad and P mo are computed from the initial conditions CR. Ta d{k) , S Ta/ ±(k)) = ((1,0), 
(0,1)), and: 

\ ~f" J iso) 

where maximally correlated (anticorrelated) modes correspond to/3 = +l(/3 = — 1). 



The extra cross-correlated term is given by 

P cm (k) = -(A;/fe ) ncor [P ad (A:)P iso (A:)]i (2.14) 

In this paper we adopt the curvaton scenario as a working example for a model that 
gives rise to a small fraction of correlated CDI isocurvature. As shown in the next sections, 
our analysis can be generalized to models with an arbitrary amount and type of isocur- 
vature. In general, by tuning the decay dynamics of a, the curvaton scenario allows for 
mixed adiabatic and isocurvature fluctuations, with any residual isocurvature perturbation 
correlated or anti-correlated to the adiabatic density one, and with the same tilt for both 
spectra. Our findings are therefore quantitative only for this case, however our conclusions 
will not be too dissimilar for a model such as the axion-like isocurvature, where rii so is fixed 
to be 1, given that n a d is not too far away from scale invariance. 

In our model we assume /3 = — 1, n s = n a d = m so and n cor = 0. 

3. Method 

In order to extend the analysis done in Ref. [39] on the impact of an isocurvature con- 
tribution on cosmological parameter estimation, we compute joint Fisher matrix forecasts 
[46] for a CMB Planck-like experiment 3 combined with a future space-based LSS survey 
with characteristics not too dissimilar from those of the proposed Euclid mission. In what 
follow we will refer to such a survey as "Euclid-like" or "LSS survey" and use the spec- 
ifications reported in [47]. We adopt the empirical redshift distribution of Ha emission 
line galaxies derived by [69] from observed Ha luminosity functions, and the bias function 
derived by [70] using a galaxy formation simulation. In particular, we choose a flux limit 
of 3xl0 -16 ergs -1 cm -1 , a survey area of 20,000 deg 2 , a redshift success rate e = 0.35, a 
redshift accuracy of cr z /(l + z) < 0.001, and a redshift range 0.45 < z < 2.05. See Tab. 1 
for a summary of the experimental specifications. 

In both cases we assume a curvaton fiducial model with mixed adiabatic and isocur- 
vature initial conditions. The fiducial value chosen for the isocurvature parameter is 
f iso = -0.08. This is well within the WMAP7-only and slightly outside the WMAP7+BAO+SN 
95%CL limits [3, 45]. On the other hand, these bounds are obtained assuming a ACDM 
model, while accounting for curvature and a varying dark energy equation of state weakens 
considerably the constraints on the isocurvature contribution [30, 31]. 

In any case it is important to stress that the Fisher analysis results, as found also by 
Ref. [39] , are not strongly affected by the choice of the fiducial values of the isocurvature 
amount within these limits. 

The set of parameters we use is: 

q a = {tt mo ,n X o,tt bo ,h, n s ,w ,w a ,log(l0 10 A tot ), f iso } (3.1) 

with fiducial values: 

q a = {0.2707, 0.7293, 0.0451, 0.703, 0.966, -0.95, 0, 1.3909, -0.08}, (3.2) 



3 http://www.sciops.esa.int/index.php?project=PLANCK 



where £l mo , Qx and fi& are, respectively, the total matter, dark energy and baryon density 
parameters at present time, h is related to the Hubble parameter by Hq = 100 h Km/(sec 
Mpc), n s and A tot are the scalar spectral index and the total amplitude accounting for 
both the adiabatic and the isocurvature contributions, wq and w a are related to the dark 
energy equation of state parameter w(z) via the standard CPL parametrization [49, 50]: 

w(z) = w + w a ——- (3.3) 

r -j- z 

Finally, fi so is the isocurvature fraction parameter of Eq. (2.5) (see §2 for detals). We 
assume flatness in the fiducial model but f2 mo and £lx are f ree parameters in the Fisher 
matrix computation: the curvature flk = 1 — ^m ~~ ^x is not fixed and it is a parameter 
in the analysis. 

Recall that, in the case of a CMB-only Planck-like experiment, the constraints on 
the isocurvature contribution are considerably weakened when f2/% is a free parameter, 
and, when also a varying dark energy equation of state is considered, the cosmological 
constraints from CMB alone weaken even more because of degeneracies [31]. 

As previously mentioned, in Ref. [39] it was shown that even a small isocurvature 
contribution, if not accounted for, can lead to an incorrect determination of the cosmological 
parameters and can affect the r s (zd) calibration from the CMB; this affects the BAO 
location estimation and, therefore, leads to systematic shifts in the derived cosmological 
parameters, shifts that result to be larger than statistical errors from forthcoming LSS 
surveys. One of the main purpose of this work is to see if this is still the case for a 
joint analysis of CMB and LSS when the information from the P(/c)-shape are included. 
Therefore, here we extend the work of Ref. [39] not only by including curvature and dark 
energy parameters in the analysis, but also by exploiting BAO location measurements 
combined to the galaxy power spectrum shape. The additional cosmological information in 
the power spectrum should help break degeneracies and thus indirectly improve constraints 
on the isocurvature parameter. 

3.1 Fisher matrix approach 

We combine the cosmological information of a CMB experiment like Planck (enclosed in 
the pCMB-Pianck Fi s h er matrix) and a galaxy survey like Euclid (enclosed in the F LSS 
Fisher matrix), so our full Fisher matrix is: 

77 _ jpCMB-Planck . r^LSS (o A \ 

*ij = i'ij + *ij ■ (3.4) 

The Fisher matrix is defined as the second derivative of the natural logarithm of the 
likelihood surface about the maximum. In the approximation that the posterior distribution 
for the parameters is a multivariate Gaussian 4 with mean /x = (x) and covariance matrix 



4 In practice, it can happen that the choice of parametrization makes the posterior distribution slightly 
non-Gaussian. However, for the parametrization chosen here, the error introduced by assuming Gaussianity 
in the posterior distribution can be considered as reasonably small, and therefore the Fisher matrix approach 
still holds as an excellent approximation for parameter forecasts. 
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100 
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217 
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Temperature sensitivity ot (^) 


2.5 
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4.8 


Polarization sensitivity o~p {^-) 
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4.2 


9.8 


Zmin 


Zmax 


Survey Area (deg 2 ) 


L^Zbin 


0.45 


2.05 


20.000 


0.1 



Table 1: TOP: experimental specifications for the Planck experiment. BOTTOM: LSS survey 
settings, i.e. redshift range z m i n and z max , survey area (A) in squared degrees, and width of the 
redshift bin Azi,i n . 



(xx*) — /x/x*, its elements are given by [51, 52, 53, 46] 
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(3.5) 



where x is a N-dimensional vector representing the data set, whose components Xj are, for 
example, the fluctuations in the galaxy density relative to the mean in N disjoint cells that 
cover the three-dimensional survey volume in a fine grid. The {0j} denote the cosmological 
parameters within the assumed fiducial cosmology. 

For the CMB Fisher matrix set up see appendix A and Tab. 1, for the LSS Fisher 
matrix set up we follow closely Ref. [41] as follows. 

3.2 Large scale structure treatment 

In order to explore the cosmological parameter constraints from a given redshift survey, 
we need to specify the measurement uncertainties of the galaxy power spectrum. The 
statistical error on the measurement of the galaxy power spectrum P g (k) at a given wave- 
number bin is [48] 



P. 



2(2tt) s 



1 



where n g is the mean number density of galaxies, Purvey is the comoving survey volume of 
the galaxy survey, and [i is the cosine of the angle between k and the line-of-sight direction 
fj, = k ■ f/k. 

In general, the observed galaxy power spectrum is different from the true spectrum, and 
it can be reconstructed approximately assuming a reference cosmology (which we consider 
to be our fiducial cosmology) as (e.g. [54]) 
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(3.7) 



where 



P g (k Ie i±,k Te q,z) = b(zY 
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matter 
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(3.8) 



In Eq. (3.7), H{z) and Da(z) are the Hubble parameter and the angular diameter distance, 
respectively, and the prefactor (Da(z) 2 c{ H (z)) / (Da(z) 2 H (z) re f) encapsulates the geomet- 
rical distortions due to a reference cosmology different from the true one [54, 57]. The 
quantities evaluated in the reference cosmology are distinguished by the subscript 'ref, 
while those in the true cosmology have no subscript. k± and k\\ are the wave-numbers 
across and along the line of sight in the true cosmology, and they are related to the wave- 
numbers calculated assuming the reference cosmology by fc rc f_i_ = k±D a(z) / D A(z) re f and 
fc re f|i = k\\H(z) re f/H(z). P s hot is the unknown shot noise that remains even after the con- 
ventional shot noise of inverse number density has been subtracted and which we assume to 
be white noise [54] . Such contribution could arise from galaxy clustering bias even on large 
scales due to local bias [58] or to misestimates of the local density In Eq. (3.8), b(z) is the 
linear bias factor between galaxy and matter density distributions, and j3(z) = f g (z)/b(z) 
is the linear redshift-space distortion parameter [59]. Thus our analysis is limited to large 
enough scales where bias is scale-independent (we will return to this below). We estimate 
fg(z) by integration of the evolution equations of the linear density perturbations. The 
linear matter power spectrum P ma .ttei(k, z) in Eq. (3.7) takes the form 
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matter 
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25/i iZ m 



tot 



T 2 (k) 



G(z) 



G{z = 0) 



-fcV<7r 



(3.9) 



where G{z) is the linear growth-factor, whose fiducial value in each redshift bin is computed 
through numerical integration of the differential equations governing the growth of linear 
perturbations in the presence of dark energy [60]. The linear transfer function T(k) depends 
on matter and baryon densities (neglecting dark energy at early times), and is computed 
in each redshift bin using CAMB 5 [61] modified to include the isocurvature contribution. 
The parameter ^4 t ot refers to the amplitude accounting for both the adiabatic and the 
isocurvature contributions as described in §2.3, and n s is the spectral index which in our 

so that it can be factorized out in the Pm&tter(k,z) expression 



n a d 



ri; 



case reads n s 
in Eq. (3.9). 

Moreover, in Eq. (3.9) we have added the damping factor exp(— k 2 /j. 2 a 2 ), due to redshift 
uncertainties, where a T = (dr/dz)a z , r(z) being the comoving distance [62, 54]. 

From the above, it should be clear that, when adding LSS, we exploit the information 
from both the galaxy power spectrum shape and BAO distance indicators. 

In each redshift shell, with size Azu n = 0.1 and centred at redshift Zj, we choose the 
following set of parameters to describe the observed power spectrum -P bs(&ref_i_, &ref|h z): 



{H(zi),DA{z i ),G(z i ),l3(z i ) 1 P l shot , uj m , u b , f iso , n s , h] 



(3.10) 



where fi so is the isocurvature fraction parameter as defined in §2, ui m = il m() h 2 , and 
ujb = £lb h 2 . Finally, since G(z), b(z), and the power spectrum normalization P are 
completely degenerate, we have introduced the quantity G(zi) = (Po)°' 5 b(zi)G(zi)/G(zo) 

[63]. 



3 http://camb.info/ 



In the limit where the survey volume is much larger than the scale of any features 
in P o bs(£0, it has been shown [64] that it is possible to redefine x n to be not the density 
fluctuation in the n th spatial volume element, but the average power measured with the 
FKP method [48] in a thin shell of radius k n in Fourier space. Under these assumptions, 
the redshift survey Fisher matrix can be approximated as [52, 64] 
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where the derivatives are evaluated at the parameter values pi of the fiducial model, and 
V e fi is the effective volume of the survey: 



Ves(k,n) 



n z Pg{k,y) 
n g P g (k,fi) + l 



•'survey) \^-^^) 



where we have assumed that the comoving number density n g is constant in position. 
Due to azimuthal symmetry around the line of sight, the three-dimensional galaxy redshift 
power spectrum P y )S (k) depends only on k and /x, i.e. is reduced to two dimensions by 
symmetry [54]. 

To minimize non-linear effects, we restrict wave-numbers to the quasi-linear regime, 
so that fc max is given by requiring that the variance of matter fluctuations in a sphere of 
radius R is cr 2 (R) = 0.25 for R = 7r/(2/c max ). This gives /c max ~ O.l/i Mpc -1 at z = 
and fc max — 0.2h Mpc -1 at z = 1, well within the quasi-linear regime. In addition, we 
impose a uniform upper limit of k max < 0.2h Mpc -1 (i.e. A; max = 0.2h Mpc^ 1 at z > 1), to 
ensure that we are only considering the conservative linear regime, essentially unaffected 
by non- linear effects. In each redshift bin we fix k m \ n = 10 -4 /i/Mpc, and we have verified 
that changing the survey maximum scale /c m j n with the shell volume has almost no effect 
on the results. 

We do not include information from the amplitude G(zi) and the linear redshift space 
distortions /3(zi), so we marginalise over these parameters and also over P l shot - Moreover, 
since we are limiting our analysis to quasi-linear scales, we do not include the effect of 
the non-linear (incoherent velocities) redshift space distortions, the so-called "Fingers of 
God" (FoG). In fact, as shown in Ref. [41], on such scales the effect of marginalisation 
over incoherent velocities can be compensated by the inclusion of growth information, and 
viceversa, if growth is marginalised over, the FoG effect is in some sense already taken into 
account, due to the tight correlation between these two effects on linear scales. 

We project p = {H(zi), DA(zi),uj m ,uJb, fiso-, n s , h} into the final sets q of cosmological 
parameters described in Eq. (3.1)-(3.2) [65, 67]. In this way we adopt the so-called "full 
P(/c)-method, marginalized over growth-information" [68], and, to change from one set of 
parameters to another, we use [65] 
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where F^§ is the survey Fisher matrix for the set of parameters q, and Fjj is the survey 
Fisher matrix for the set of equivalent parameters p. 

We combine the LSS Euclid-like survey with CMB information from the Planck satel- 
lite, using the parameter set as specified in Appendix A. 

After marginalization over the optical depth, we propagate the Planck CMB Fisher 
matrix F i ■ ~ anc into the final sets of parameters q, by using the appropriate Jacobian 
for the involved parameter transformation. 

The 1-<t error on q a marginalized over the other parameters is <j(q a ) = \/(F~ 1 ) aa , 
where -F -1 is the inverse of the Fisher matrix. We then consider constraints in a two- 
parameter subspace, marginalizing over the remaining parameters, in order to study the 
covariance between fi so and other relevant parameters. 

3.3 Correlation and shifts 

To quantify the level of degeneracy between the different parameters, we estimate the 
so-called correlation coefficients, given by 



(.-^ )PaP/3 



{■T )PaPa K* 1 >PPPP 



(3.14) 



where p a denotes one of the model parameters. When \r\ = 1, the two parameters are 
totally degenerate, while r = means they are uncorrelated. 

Furthermore, we use another powerful tool encoded in the Fisher analysis which allows 
to calculate the shift in the best fit parameters 59 a in the case that the isocurvature amount 
parameter fi so is fixed to a wrong fiducial value, without recomputing the full covariance 
matrix. 

In general, given a number p of parameters \I/ 7 (7 = 1, ..,p) fixed to an incorrect value, 
which differs from the "true" value by an amount 5^ 7 , the resulting shift on the best fit 
value of the other n parameters 9 a (a = 1, ..., n) is (e.g. [66]): 

89 a = -{F- 1 ) a$ Sfr6V 1 . (3.15) 

Here (-F -1 ) B is the sub-matrix of the inverse Fisher matrix in Eq. (3.4), correspond- 
ing to the 9 a parameters (i.e. the inverse of the Fisher matrix, without the rows and 
columns corresponding to the "incorrect" parameters), and Sg 7 is the Fisher sub-matrix 
including also the \I/ 7 parameters. The shifts in the best fit parameters induced by 
setting the amount of isocurvature fi so to an incorrect value is: <5\I/ 7 = 5fi so where 
9 a = {ttm ,Qx ,tt bo ,h,n s ,w ,w a ,log(W 10 A tot )}. 

4. Results 

In this work we consider a cosmological model with curvature, isocurvature and varying 
dark energy. In this case, the combination of a Planck-like CMB experiment and a LSS 
probe improves the constraints on the cosmological parameters by at least one order of mag- 
nitude with respect to the Planck-only constraints for a less general model with curvature 
and isocurvature but with a cosmological constant (see e.g. [36]). 
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Figure 1: 2-parameter confidence levels for fi SO and q a with q a = h, fif,, f2 m ,log(10 10 A tot ) for 
the fiducial model with f iso = —0.08 obtained for an Euclid-like survey only (larger contours) and 
after combining the survey data with Planck. For log(10 10 A to t) we show the constraints from 
LSS+CMB only since for the LSS alone, using the P(£;)-method, we marginalize over the power 
spectrum normalization. The blue dotted line and the red dashed line represent the 68% C.L., 
95.4% C.L., respectively. The black solid line shows the 1-parameter confidence level at 1-a. 



Here, we have extended the analysis of Ref. [39] to the more general model described 
in §3, and, moreover, we have added the cosmological information from the LSS to the 
CMB-only constraints. The results are summarized in the first three columns of Tab. 2. 
Note that the 1-a error on fi SO obtained from the joint CMB+LSS analysis is ~ 0.008; this 
is slightly better than the error expected from Planck alone, <t(/j so ) = 0.01 of Ref. [39], 
where 0^ = and w = — 1 (both kept as fixed parameters) were assumed. 

In Figs. 1-2 we show the contours at 68% CL, 95% CL and 1-parameter CL at 1-<t, 
for fi so and q a with q a = {n s , Qxo, wo,w a , h, O^, Q m , log(10 10 A tot )}, for the fiducial model 
with fi so = —0.08, obtained from an Euclid-like survey alone (larger contours), and from 
the combination with Planck data. For log(10 10 A tot ) we show only the constraints from 
LSS+CMB, since for the LSS alone, using the P(/c)-method, we marginalize over the power 
spectrum normalization. 

From Figs. 1-2 and column 3 of Tab. 2 it is evident that the CMB+LSS constraints 
improve dramatically with respect to the LSS-only case and also with respect to the CMB- 
only case [39, 36]. This means that adding LSS is a powerful tool which allows to break 
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Figure 2: 2-parameter confidence levels for fi SO and q a with q a = n s ,flxo,wo,w a for the fiducial 
model with fi SO = —0.08 obtained for an Euclid-like survey only (larger contours) and after com- 
bining the survey data with Planck. The blue dotted line and the red dashed line represent the 
68% C.L., 95.4% C.L., respectively. The black solid line shows the 1-parameter confidence level at 
1-cr. 



the strong degeneracies that affect CMB measurements when curvature, isocurvature and 
varying dark energy are included in the analysis. As a consequence, constraints on the 
nature of the initial conditions are greatly improved by the combination CMB+LSS, with 
respect to CMB alone. 

The next issue we address regards possible systematic errors on cosmological parame- 
ters due to incorrect assumptions about the nature of the initial conditions. 

To this purpose, we calculate the shifts on the cosmological parameters due to an 
incorrect amount of the isocurvature fraction fi so in the fiducial model. 

The results are summarized in the last three columns of Tab. 2. As can be noticed, 
the absolute shifts (see Eq. (3.15)) 59 a /5fi so assume the largest values for the dark energy 
parameters Ox , Wq, w a , for the spectral index n s , and for the amplitude A tot (see the 4 th 
column in Tab. 2). 

In interpreting these results one must bear in mind that, rather than the shifts theme- 
selves, it is important to determine if and in which cases these shifts on cosmological 
parameters are larger than the corresponding statistical errors (3 column in Tab. 2). 

We have assumed a fiducial value fi so = —0.08, so, if the true underlying model were 
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Table 2: Forecast 1— a errors of the cosmological parameters q a considered in the text and cor- 
responding correlations with fi SO , for an Euclid-like experiment alone and in combination with 
Planck. In particular, in the 2 nd column we report the fiducial values adopted for the parameters 
listed in the 1 st column. The corresponding \-a errors are reported in the 3 rd column, and the 
q a ~fi SO correlation coefficients in the 4 th column. Finally, in the 5 th , 6 th and 7 th columns we list 
the shifts of the parameters 9 a for Sf iso — 1,0.08,0.01, respectively. 



adiabatic (for which the "true" value for the isocurvature fraction is fi so = 0), fi so would 
be fixed to an incorrect value by 5fi so = +0.08. In this case, the correspondent shifts on 
the cosmological parameters are shown in the 5 column of Tab. 2. For O, mo , £lx , ^tot 
and n s (this last being the most affected parameter) the shifts are bigger than the 1—cr 
statistical errors, meaning that the corresponding degeneracies with fi so are not completly 
solved, and a wrong assumption on the amount of isocurvature perturbations produces a 
bias in the estimation of the cosmological parameters even when adding LSS to CMB. On 
the other hand, the shift is slightly smaller but still comparable to the 1—a error for w a . 
For the case 5fi so = +0.08, adding the LSS information breaks the degeneracy fi so -h with 
respect to the CMB-only case, even if the /i-shift is still comparable to the statistical error. 
For the other parameters the degeneracies are reduced and the shifts are smaller but still 
of the same order of magnitude as the statistical \-a errors. 

On the other hand, if we consider the smaller value 5fi so = +0.01, the LSS informa- 
tion from an Euclid-like survey solves the degeneracies with fi so for all the cosmological 
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parameters except n s , as shown in the 6 column of Tab. 2. This is a remarkable result 
given that the analysis is done including curvature and a varying dark energy equation of 
state. 

Since the statistical errors are not very sensitive to the fiducial fi so value, the case 
5 fiso = +0.01 is representative to quantify the difference in modeling with the pure adia- 
batic case and the mixed underlying model with a fiducial fi so = —0.01. 

As shown in Figs. 1-2, for CMB+LSS the constraints on parameters which are degen- 
erate with fi so improve considerably. However, note that the shifts on parameters as £lx, 
fi m o, wq and w a are larger than in the LSS-only case, as shown in column 5 of Tab. 2. 

This is confirmed by the behavior of the correlation coefficient r (see Eq. (3.14)), shown 
in column 4 of Tab. 2. For instance, if we consider the correlation fi so -^x , we can observe 
that fi so and £lx are negligibly correlated for LSS measurements, but they are tightly 
correlated for CMB measurements; when combining the two information, the LSS help to 
alleviate the degeneracy but a level of correlation still remains. 

In some cases the correlation coefficient r, as well as the shifts, can even change sign 
when adding Planck to LSS. This change follows the dominant correlations between fi so 
and the cosmological parameters either from LSS or CMB. 

The parameter for which the degeneracy with the isocurvature fraction parameter /j so 
is not solved even when combining the CMB and LSS information is the spectral index 
n s . We interpret this as due in part to the fact that in this paper we assume a curvaton 
underlying model as a working example for which n s = n a( i = rii so . Therefore, n s encodes 
the information of both the spectral indices which are correlated with fi so for both CMB 
and LSS. In fact, the effect of a variation of rii so on the C\ and the P(k) would be different 
even for a fixed fi so - Conversely, the effect of changing fi so would be different in CMB 
and LSS if n ac i were not to be equal to rii so - For the current data, the effect of changing 
fi so is stronger in CMB which thus "dominates" over LSS. Therefore, the magnitude of 
the degeneracy n s — fi so is related to our choice n s = n a( i = rii so - When rii so = n a d the 
information is compressed in one parameter which in some way is "double-degenerate" with 
fiso and adding LSS helps but not enough to fully solve the degeneracy. 

These results confirm once again the great potential of forthcoming large scale galaxy 
surveys. 

5. Discussion and conclusions 

In this work we have relied on the Fisher matrix approach to forecasts, which is not well 
suited to quantify systematic effects. The errors we have reported are statistical only, which 
are optimistic in the sense that we assume that systematic errors are vastly sub-dominant. 
Let us therefore discuss sources of systematic errors of major concern and their possible 
consequences on our results: these are the effects of non-linearities and galaxy bias. 

In the analysis presented here we have used the linear theory matter power spectrum 
and a scale-independent galaxy bias. Our choice of the maximum k-mode ensures that the 
effects of non-linearity are quite mild and that eventually, before an analysis is performed 
on real data, they could be accurately modeled. Having considered only the mildly non- 
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linear regime ensures that non-linearities will not erase or reduce significantly the signal. 
Reconstruction techniques that are being developed [71, 72] and/or are already being ap- 
plied to surveys [73] offer promising prospects. Thus while it will be mandatory to include 
non-linearities in the actual data analysis, the forecasted errors are not made artificially 
smaller by using the linear matter power spectrum to compute our Fisher matrices. 

Including smaller, more non-linear scales might help to improve the errors but would 
worsen the systematic effect of non-linearities and of a possible scale-dependent and/or 
non-linear bias. In this work we have made the simplifying assumption that bias is scale- 
independent but that the redshift dependence is not well known and so it is marginalised 
over. Bias is known to be scale-independent on large scales but scale dependent for small 
scales and/or very massive halos. While halo bias can be modeled, the bias of galaxies 
hosted in dark matter halos can be more complicated especially at small scales. Here we 
have considered scales where the so-called two halo term dominates and so complicated 
bias arising from the details of galaxy formation and halo occupation distributions do 
not matter. In addition we have checked that most of the information, useful to break 
degeneracies between the sound horizon at radiation drag and isocurvature contribution, 
come mostly from very large scales (those constraining the matter -radiation equality) . The 
BAO signal is localized on scales where scale-dependent bias is more likely to be an issue, 
but the BAO location is relatively robust to this. Our reported error-bars could have been 
made smaller by considering also information about the growth of structures (i.e., redshift 
space distortions and the normalization of the power spectrum), however we decided to, 
conservatively, marginalise over the amplitude to be less sensitive to possible assumptions 
about galaxy bias. 

As found in Ref. [39] for a Planck-like CMB experiment, even a tiny isocurvature 
fraction contribution, if not accounted for, can lead to an incorrect determination of the 
cosmological parameters even for a standard flat model with a cosmological constant. In 
this case however, when the isocurvature parameter fi so is added to the modelling, the 
errors increase but the cosmological parameters are recovered correctly. The error degra- 
dation from a CMB-only analysis is particularly severe for models where spatial flatness is 
not imposed and curvature is left as a free parameter. 

In this work we extend the analysis done in Ref. [39] and present the results from a 
joint analysis of a CMB Planck-like experiment and a LSS EULID-like galaxy survey, for a 
cosmological model that includes an isocurvature fraction in the initial conditions, varying 
dark energy and curvature. 

As shown in Figs. 1-2 and in Tab. 2, we find that adding the cosmological information 
coming from an Euclid-like LSS probe strongly reduces the degeneracies introduced by the 
presence of an isocurvature contribution with respect to a Planck-like CMB probe alone. 

When isocurvature is included in the analysis, we find that the most affected parameter 
is the spectral index n s for which the degeneracy with the isocurvature fraction parameter 
fi so is never fully solved even when combining CMB and LSS. This is because in this paper 
we assume a curvaton underlying model as a working example for which n s = n a ^ = rii so . 
Therefore n s encodes the information of both the spectral indices which are correlated 
with fi so for both CMB and LSS. In this case, from the combination CMB+LSS, we 
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find an absolute n s -shift= 0.277 which is larger than the statistical 1— a error by two 
order of magnitude. For all the cosmological parameters the correlations with /j so and 
the corresponding shifts are reported in Tab. 2. Our findings are quantitative only for the 
n ad = n iso case, but can hold qualitatively for other models with rii so ~ 1, as for example 
the axion-like isocurvature model. In this case, since n ac i is close to scale invariant, we 
expect that the effect will be quantitatively very similar. 

We investigate also the effect on the recovered cosmological parameters in the case 
of the presence of a small amount of isocurvature which is however ignored in the anal- 
ysis. Ref. [39] found that even a tiny isocurvature fraction, if neglected, could bias the 
determination of "standard rulers" like the sound horizon r s {z ( i) at radiation drag and 
introduce systematic biases in the interpretation of BAO observables that are larger than 
the statistical error-bars forecast for future and forthcoming surveys. Ref. [39] considered 
a flat Universe where the dark energy is a cosmological constant, but in a more general 
cosmology the effect can be dramatically worst. 

Here we find that adding LSS information to CMB priors eliminates the bias (i.e. 
reduces it well below the statistical errors) in the case of a small amount of isocurvature 
fraction fi so ~ 0(1O~ 2 ), even for a general cosmological model with varying dark energy 
and curvature. The only parameter for which the systematic shift remains still comparable 
to the 1— a error is the spectral index n s . 

These remarkable results are obtained by better exploiting the potential cosmological 
information of future LSS surveys by adding to BAO measurements also the power spec- 
trum shape, following the so-called P(/c)-method marginalised over growth information. In 
fact, the additional cosmological information in the power spectrum helps to break degen- 
eracies among the cosmological parameters, and thus indirectly improves the constraints 
on the isocurvature parameter fi so . 

Considering that -conservatively- we have not included information on the power spec- 
trum amplitude and on the growth of cosmological structures (we marginalize over redshift 
space distortions parameters and LSS power spectrum normalization), we conclude high- 
lighting once again the synergy and complementarity between CMB and LSS and the great 
potential of future galaxy surveys. 
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A. Planck Fisher matrix 

In this work we use the Planck mission parameter constraints as CMB priors, by estimating 
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the cosmological parameter errors via measurements of the temperature and polarization 
power spectra. As CMB anisotropics, with the exception of the integrated Sachs- Wolfe 
effect, are not able to constrain the equation of state of dark energy (wo,w a ) 6 , we follow 
the prescription laid out by DETF [56]. 

We do not include any B-mode in our forecasts and assume no tensor mode contribu- 
tion to the power spectra. We use the 100 GHz, 143 GHz, and 217 GHz channels as science 
channels. These channels have a beam of 6>f w h m = 9.5', 0f w h m = 7.1', and 0f w h m = 5', 
respectively, and sensitivities of gt = 2.5fiK/K, ot = 2.2fxK/K, gt = A.8fiK/K for 
temperature, and up = 4fiK/K, op = 4.2fiK/K, op = 9.8fj,K/K for polarization, respec- 
tively. We take / s k y = 0.80 as the sky fraction in order to account for galactic obstruction, 
and use a minimum ^-rnode £ m ; n = 30 in order to avoid problems with polarization fore- 
grounds and not to include information from the late Integrated Sachs- Wolfe effect, which 
depends on the specific dark energy model. We discard temperature and polarization data 
at £ > 2000 to reduce sensitivity to contributions from patchy reionization and point source 
contamination (see [56] and references therein). 

We assume a fiducial cosmology with an anti-correlated isocurvature contribution, 
varying dark energy and curvature. Therefore we choose the following set of parame- 
ters to describe the temperature and polarization power spectra 6 = (w m , uib-, fiso, 100 x 
0s, log(10 10 A tot ) , ns , wo, w a ), where 9s is the angular size of the sound horizon at last scat- 
tering and wq and w a are the dark energy parameters according to the CPL parametrization 
of the dark energy equation of state w{z) = wo + w a z/(l + z). 

The Fisher matrix for CMB power spectrum is given by [74, 75]: 

r?CMB-Planck ST^ V^ 9Cx,l nrw -l 9Cy,l , . > 

^ =Z^^- C0V xy^-, (A.l) 

I X,Y l 3 

where 6i are the parameters to constrain, Cx,i is the harmonic power spectrum for the 
temperature-temperature {X = TT), temperature- E-polarization (X = TE) and the E- 
polarization-E-polarization {X = EE) power spectrum. The covariance COV^ Y °f the 
errors for the various power spectra is given by the fourth moment of the distribution, 
which under Gaussian assumptions is entirely given in terms of the Cxi with 

COV T , T = ft {C T ,i + W- l BY 2 f (A.2) 

COX E ^E = ft (C E ,i + Wp'Bf 2 ) 2 (A.3) 

COX T e,te = ft [C$ E ,1 + (Ct,i + W^Bf 2 ) {C E ,i + W^Bf 2 ) ] (A.4) 

COX T: e = ftC^ E ,i (A.5) 

COV TjT £ = f e C TE ,i (C T) i + W T l B- 2 ) (A.6) 

CON e ,te = feC TE ,i {C E ,i + Wp l B{ 2 ) , (A.7) 

where f t = 2/ \{2l + l)/ sky ), W t ,p = E c W t,p> W t,p = Kp^fwhm)" 2 bein S the wei S ht 
per solid angle for temperature and polarization respectively, with a 1— a sensitivity per 



6 On the contrary, using (wo,Wa) as model parameters to compute the CMB Fisher matrix could artifi- 
cially break exiting degeneracies. 
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pixel of Oj>p an d a beam of 0^ whm extent, for each frequency channel c. The beam window 
function is given in terms of the full width half maximum (fwhm) beam width by Bj = 
E c m 2w T,p/ W T,P, where (Btf = exp (-£(£ + l)/(/ s c ) 2 ), l c s = (^wh m ) _1 ^(81n2) and 
/sky is the sky fraction [76]. 

We then calculate the Planck CMB Fisher matrix with the help of the publicly available 
CAMB code [77]. Finally, we transform the Planck Fisher matrix for the DETF parameter 
set to the final parameter sets q considered in this work (see §3), using the transformation 



j?C MB- Planck _ \~^ uu i rpCMB-Planck uu j i a c 
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